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The flowfield characteristics in a rocket engine coolant channels are analyzed by use of a 
commercial CFD and multiphysics software developed by the CFD Research Corp. called 
CFD-ACE+. The channels are characterized by high Reynolds number flow, varying aspect 
ratio, varying curvature, asymmetric and symmetric heating. The supercritical hydrogen 
coolant introduces large property variations that have a strong influence on the developing 
flow and the resulting heat transfer. This paper only shows the effect of aspect ratio and 
curvature for constant properties. 
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Nomenclature 

Pressure Drop 

Bulk Inlet Temperature 

Bulk Outlet Temperature 

Outlet Wall Temperature 

Density 

Viscosity 

Thermal Conductivity 
Constant Pressure Specific Heat 
cylinder diameter 
Perimeter 
Channel Length 
Mass Flow Rate 
Velocity 

Total Heat Transfer Rate 
Constant Surface Heat Flux 

Heat Transfer Coefficient 
Friction Factor 
Reynolds Number 
Prandtl Number 
Nusselt Number 


I. Introduction 

Because of high combustion temperatures and high heat transfer rates from the hot gases to the walls of the 
thrust chamber, cooling is a major design consideration in a rocket engine. There are many types of cooling. One of 
the most common and the one of particular interest in this paper is regenerative cooling. Regenerative cooling is 
normally used to protect the walls from the severe thermal environment that it sees during the combustion process. 
In most regeneratively cooled engines, the coolant passages are rectangular in shape. These tubes run the length on 
the combustor with the coolant entering the supersonic end, passing through the throat of the combustor, and exiting 
near the injector face. Liquid hydrogen is a typical coolant, and the pressures seen inside the thrust chamber are 
sufficient to ensure that the fluid remains supercritical throughout the combustion process. 
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Even though the designs of these cooling passages are critical for success, their heat transfer characteristics are 
not fully understood. The heat transfer characteristics of rectangular tubes differ in many ways to those of circular 
tubes. One clear difference is the way the fluid is heated. Heat transfer research often entails the use of an 
electrically heated circular tube, which is not representative of typical applications. Most regen tubes are heated 
from one side only. The combination of highly conductive materials and the milled rectangular passages that are 
used in design of the combustor body provides a fin effect between adjacent coolant passages which makes the 
rectangular cavities more effective. Also the curvature introduced by the shape of the combustor generates vortices 
inside the channels which enhance the heat transfer. Another factor is the strong compressibility of supercritical 
hydrogen. This strong compressibility will certainly have an effect on the flowfield. With all these different 
phenomena that occur inside the tube, the aspect ratio that is used in the combustor design represents a key factor 
that can be used to optimize the transfer of heat from the chamber walls to the coolant fluid. 

To understand the heat transfer environment inside these rectangular tubes, the simple case of a straight 
rectangular tube is considered. Even in straight channels the flow field remains complex due to the fin effect, high 
Reynolds number, and the supercritical coolant. The supercritical conditions inside the coolant passages of rocket 
engines yield fluid properties that are strong functions of both pressure and temperature. The specific property 
variations that must be considered are specific heat, thermal conductivity, density and viscosity. The current study 
will only deal with constant properties. Variable property comparisons are left to future investigations. 

II. Code Validation 

In any numerical simulation, questions directed to the accuracy of the results always arise. Because of this, code 
validation is a must to give credibility to the results. Since we are dealing with a high Reynolds number flow, the 
turbulent y+ value must be taken into account along with the grid mesh to insure that proper flow characteristics are 
being simulated. Also, since we are heating the tubes, validation of the temperature prediction must be made. 

Although we are studying the heat transfer effects in rectangular passages with high aspect ratio, with and 
without curvature, we first need validate the code using a simple case that is well known. One case that is well 
documented is circular tubes with constant surface heat flux. Temperature and pressure drop predictions are very 
important for the design of regenerative cooling passages. Therefore, if the commercial software can predict the 
pressure drop, bulk outlet temperature and surface temperature with some accuracy, then the results obtained for 
rectangular passages can be accepted with some degree of comfort. Typical tube conditions for regenerative cooled 
engines were used. The tube diameter used is 2 mm with a tube length of 34.25 mm. The inlet fluid temperature is 
set to 40 K, with an inlet mass flow rate of 5.516E-06 kg/s. A constant surface heat flux 100,000 W/m 2 * is applied to 
the wall. Thermal physical properties for supercritical hydrogen at 40 K, such that the Reynolds number is 10,000 is 
taken from Ref. 1. To calculate the outlet temperature, the following energy balance equation is employed along 
with the constant surface heat flux equation 2 . 


(Icon, =™ C p (T mo -T mi ) 


< 1co m = ^*L) 

Utilizing equations (1) and (2), the bulk outlet temperature equation is 

T 

m.o . m,i 

mc p 

Now to find the surface temperature at the outlet, we need Newton’s Law of cooling 2 . 

q>KT„ -T m .) 

To find /i, the turbulent Nusselt number correlation by Dittus-Boelter is employed 2 . 
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To find the pressure drop, we need the following equation derived from the Darcy-Weisbach equation and the 
friction factor/ from the Moody chart 3 . 


A P = f 


pLV 2 

2D 


( 6 ) 


From the Moody chart, /is approximately 0.0309 for a smooth tube of Reynolds number 10,000. The results of 
the theoretical and numerical calculation are shown in Table 1 below. 



Bulk Outlet Temperature 

(K) 

Outlet Wall Temperature 

(K) 

Pressure Drop 
(Pa) 

Theoretical 

51.08 

116.43 

2.899 

CFD 

53.82 

112.7 

3.208 


Table 1: Numerical and Theoretical comparison 


For our numerical simulation, we used a k-s turbulence model, a turbulent Prandtl number of 0.9, and a fixed 
pressure outlet boundary condition. The K-e turbulence model is a standard turbulence model. It is a high Reynolds 
number model that uses wall functions for cells that are adjacent to the wall. CFD Research Corp. recommends a 
turbulent y+ value between 30 and 150 for the first cell adjacent to the wall. The mesh used is shown if Figure 1 and 
Figure 2. The grid that is shown in Figure 1 is a cross section of the pipe and shows grid stretching away from the 
wall. Figure 2 shows the grid for the length of the pipe and is stretched toward the inlet to aid in flow development. 
This mesh consists of 50,000 nodes. Also, we set the turbulent kinetic energy according to the turbulent intensity of 
0.051 and a specific dissipation rate based on hydraulic diameter 0.002. 

In Figure 3, we can see the streamwise contours of pressure, velocity and temperature. From Table 1, we can see 
that the numerical model prediction is close to that of the analytical solution. The numerical model varies less than 
6% for the bulk outlet temperature and the outlet surface temperature than that of the analytical solution. The 
pressure drop for the code is around 10.7% off from the analytical solution. Overall, these differences are small and 
we can feel confident about the numerical model. 



Figure 1: Computational grid in the cross plane for Reynolds number 10,000. 
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Figure 3b 
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Figure 3c 

Figure 3: Streamwise plots of the pressure, velocity, and temperature down the center plane. 


III. Results and Future Plans 

Now that the code has been show to be capable of predicting turbulent channel flow, we can now proceed to 
looking at several different cross-sectional shapes and channel curvature. This parametric study of aspect ratio holds 
the cross sectional area fixed. Also, the length of the tube considered will remain fixed. Table 2 shows a list of 
different aspect ratios and curvatures and the current status of each trial. 



Reynolds number 10,000 

Aspect Ratio 

Straight 

600/300 

90 Deg bend 
600/300 

Straight 

300/300 

90 Deg bend 
300/300 

1.5 to 1 

Complete 

Complete 

Complete 

Complete 

5 to 1 





8 to 1 





10 to 1 






Reynolds number 500,000 

Aspect Ratio 

Straight 

600/300 

90 Deg bend 
600/300 

Straight 

300/300 

90 Deg bend 
300/300 

1.5 to 1 





5 to 1 





8 to 1 





10 to 1 






Table 2: Trial status table 


Since the overall progress is somewhat in the beginning stages, not much can be said about how aspect ratio 
effects heat transfer and pressure drop. Figure 4 and Figure 5 are some temperature plots near the exit from the 
straight and 90 degree bend channels with wall temperature of 600 K and 300K. 
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Figure 5: 90 Bend Channel Near Exit Plane Temperature Contour 
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IV. Conclusion 

It is well known that there exists a secondary flow in high aspect ratio tubes with curvature 6 . With this secondary 
flow, more heat can be transferred from the hot combustor wall to the coolant channel. From Figure 4 and Figure 5 
you can start to see some of the effects of this secondary flow. We have cooler liquid reaching the hot combustor 
wall and adding more heat into the centerplane shown in Figure 5. It is expected that as the aspect ratio increases, 
the more effective the cooling process will become. 


7 

American Institute of Aeronautics and Astronautics 


Acknowledgments 

The work reported in this paper was supported by a contract from the National Aeronautics and Space 
Administration through the Constellation University Institutes Project (CUIP) and the Rocket Engine Advancement 
Program (REAP). 


References 

‘Linstrom, P. J. and Mallard, W.G., Eds., NIST Chemistry WebBook, NIST Standard Reference Database Number 
69, National Institute of Standards and Technology, Gaithersburg, MD, June 2005. ( http://webbook.nist.gov ). 
2 Incorpera, Frank P., and DeWitt, David P., Fundamentals of Heat and Mass Transfer , 5 th ed., John Wiley &Sons, 
New York, 2002, Chap 8. 

3 White, Frank M., Fluid Mechanics , 4 th ed., McGraw-Hil, New York, 1999, Chap 6. 

4 Yagley, J., Feng, J., and Merkle, C. L., “The effect of Aspect Ratio on the Effectiveness of Combustor Coolant 
Passages,” AIAA Paper 92-3153, AIAA/SAE/ASME/ASEE 28 th Joint Propulsion Conference and Exhibit, 
Nashville, TN .,July 6-8, 1992. 

5 Yagley, J., Feng, J., and Merkle, C. L., “CFD Analyses of Coolant Channel Flowfields,” AIAA Paper 93-1830, 
AIAA/SAE/ASME/ASEE 29 th Joint Propulsion Conference and Exhibit, Monterey, CA., June 28-30, 1993. 

6 Meyer, Michael L., “Visualization of Secondary Flow Development in Huge Aspect Ratio Channels with 
Curvature,” AIAA Paper 94-2979, AIAA/SAE/ASME/ASEE 30 th Joint Propulsion Conference, Indianapolis, IN., 
June 27-29, 1994. 

7 Meyer, Michael L., “Electrically Heated Tube Investigation of Cooling Channel Geometry Effects,” AIAA Paper 
95-2500, AIAA/SAE/ASME/ASEE 31 th Joint Propulsion Conference and Exhibit, San Diego, CA., July 10-12, 
1995. 

8 Humphrey, J. A. C., Taylor, A. M. K., and Whitelaw, J. H., “Laminar Flow in a Square Duct of Strong Curvature,” 
Journal of Fluid Mechanics, Vol. 83, Part 3, 1977, pp. 509-527. 


8 

American Institute of Aeronautics and Astronautics 


